AO- A063  723  PURDUE  UNIV  LAFAYETTE  IND  PROJECT  SQUID  HEADQUARTERS  F/G  21/2 

SECOND-ORDER  CLOSURE  MODEL IN6  OF  VARIABLE  DENSITY  TURBULENT  FLO— ETC<U) 
MAR  79  A  K  VARMA *  P  J  MANSFIELD*  G  SANDRI  N00014-75-C-1143 

JNCLASSIF1ED  SQUID-ARAP-2-PU  NL 


.  .  . 

MICROCOPY  RESOLUTION  TEST  CHART 
NATiONAL  BUREAU  OE  STANDARDS -1963-^. 


ADA083723 


PROJECT  SQUID 

TECHNICAL  REPORT  ARAP-2-PU 

SECOND-ORDER  CLOSURE  MODELING  OF 
VARIABLE  DENSITY  TURBULENT  FLOWS 


ASHOK  K.  VARMA,  PETER  J.  MANSFIELD  and  GUIDO  SANDRI 
AERONAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON,  INC. 
50  WASHINGTON  ROAD,  P.O.  BOX  2229 
PRINCETON,  NEW  JERSEY  08540 


MARCH  1979 


Projact  SQUID  if  a  cooperative  program  of  basic  research  relating  to  Jet  Propul¬ 
sion.  It  is  sponsored  by  the  Office  of  Naval  Research  and  is  assisted  by  Purdue 
University  through  Contract  N00014-75-C-1143,  NR -096- 038. 


QJ 


UJ. 


Published  for  ONR  by 
School  of  Mechanical  Engineering 
Chaffee  Hall 
Purdue  University 
West  Lafayette,  Indiana  47907 


I 

s 


DTIC 

ELECTE 
APR  2  9  198 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


80  4  29  043 


\ 


Technical  Report  ARAP-2-PU 


PROJECT  SQUID 


A  COOPERATIVE  PROGRAM  OF  FUNDAMENTAL  RESEARCH 
AS  RELATED  TO  JET  PROPULSION 

OFFICE  OF  NAVAL  RESEARCH,  DEPARTMENT  OF  THE  NAVY 

*•* 

CONTRACT  N00014flC-!143,  NR-098-038 


SECOND-ORDER  CLOSURE  MODELING  OF 
VARIABLE  DENSITY  TURBULENT  FLOWS 


by 


Ashok  K.  Varma1,  Peter  J.  Mansfield2  and  Guido  Sandri3 
Aeronautical  Research  Associates  of  Princeton,  Inc. 

50  Washington  Road,  P.O.  Box  2229 
Princeton,  New  Jersey  08540 


March  1979 


Published  for  ONR  by 
School  of  Mechanical  Engineering 
Purdue  University 
West  Lafayette,  Indiana  47907 


This  document  has  been  approved  for  public  release  and  sal 
its  distribution  is  unlimited. 


Consultant 
2Senior  Staff 
3Senior  Consultant 


SUMMARY 


Mixing  and  chemical  reactions  under  turbulent  flow  conditions  are  a 
basic  feature  of  the  energy  release  processes  in  many  combustion  and  pro¬ 
pulsion  systems.  The  development  of  predictive  calculation  procedures 
for  these  systems  requires  the  understanding  and  modeling  of  coupling 
between  turbulence  and  various  physical  and  chemical  processes.  Second- 
order  closure  modeling  of  turbulent  flows  provides  a  rational  framework 
for  studying  these  interactions. 

Models  for  the  scalar  probability  density  function  (pdf)  have  to  be 
developed  to  achieve  closure  of  turbulent  transport  equations  for  mixing 
and  reacting  flows.  A  delta  function  "typical  eddy"  pdf  model  for  two 
species  flows  has  been  developed  and  incorporated  into  a  complete  second- 
order  closure  computer  program.  The  program  has  been  used  to  study  uniform 
and  variable  density  flowfields  and  the  model  predictions  have  been  com¬ 
pared  to  experimental  measurements.  The  modeling  of  turbulence  dynamics 
for  variable  density  flows  requires  further  improvement.  However,  the 
importance  of  modeling  the  higher-order  scalar  correlations  has  been  demon¬ 
strated. 

A  number  of  statistical  constraints  on  three  species  flowfields  have 
also  been  derived.  These  will  be  useful  in  the  development  of  the  "typical 
eddy"  pdf  model  for  reacting  flows. 
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1.  INTRODUCTION 


Mixing  and  chemical  reactions  under  turbulent  flow  conditions  are  a 
basic  feature  of  the  energy  release  processes  in  gas  turbine  combustors 
and  other  propulsion  systems,  and  any  predictive  modeling  of  the  flowfield 
has  to  properly  account  for  the  effects  of  the  turbulence  of  the  flow  and 
its  interactions  with  various  physical  and  chemical  processes.  There  are 
significant  interaction  effects  of  the  turbulence  on  combustion  and  of 
the  combustion  on  the  turbulence,  which  are  important  in  determining  com¬ 
bustion  efficiency,  pollutant  formation,  combustion  noise,  heat  transfer, 
etc.  Second-order  closure  modeling  of  turbulent  flows  provides  a  conve¬ 
nient  framework  for  studying  these  interactions  and  holds  promise  of  pro¬ 
viding  a  reliable  predictive  computational  tool  for  the  design  of  new 
systems  and  improvement  of  existing  combustion  systems. 

The  modeling  of  turbulent  flows  with  variable  density  and  combustion 
requires  the  development  of  models  for  third-  and  higher-order  scalar 
correlations  or,  alternately,  a  model  for  the  joint  scalar  probability 
density  function  (pdf).  The  development  of  a  satisfactory  scalar  pdf 
model  is  one  of  the  major  objectives  of  this  research  effort.  A  "typical 
eddy"  pdf  model  composed  of  delta  functions  has  been  proposed,  and  a  model 
for  two  species  variable  density  mixing  flows  has  been  constructed  and 
extensively  tested  by  comparison  with  measured  pdf's  in  a  variable  density 
shear  layer  flow  (Varma  et  al.,  1978). 


This  pdf  model  has  been  incorporated  in  a  complete  second-order  clo¬ 
sure  computer  program  developed  at  A.R.A.P.  The  code  is  designed  for 
multispecies  turbulent  reacting  flows  and  solves  a  set  of  27  coupled, 
parabolic  partial  differential  equations  for  the  mean  variables  and  vari¬ 
ous  independent  second-order  correlations.  The  second-order  closure 
program  with  the  scalar  pdf  model  has  been  used  to  study  uniform  and 
variable  density  planar  shear  layer  flows  and  a  variable  density  axisym- 
metric  jet  flowfield.  The  program  predictions  are  compared  to  experi¬ 
mental  measurements  in  this  report  (Section  2),  and  the  importance  of 
modeling  the  higher-order  scalar  correlations  has  been  demonstrated. 

The  development  and  testing  of  a  delta  function  pdf  model  for  three 
species  is  also  continuing.  Various  statistical  constraints  on  corre¬ 
lations  for  three  species  that  have  been  derived  up  to  now  are  presented 
in  Section  3  of  this  report. 


2.  MODEL  VERIFICATION 


Uniform  and  Variable  Density  Planar  Shear  Layers 

Detailed  measurements  of  the  mean  velocity  u  ,  mean  density  p  , 
mean  species  mole  fraction  c  ,  species  correlation  c'c1  ,  and  species 
probability  density  functions  (pdf's)  in  uniform  density  and  variable  den¬ 
sity  planar  shear  layers  have  been  made  by  Prof.  Roshko  and  his  colleagues 
at  Cal  Tech.  The  data  has  been  reported  by  Brown  and  Roshko  (1974)  and 
Konrad  (1976).  The  data  appears  to  be  the  only  comprehensive  set  of 
measurements  for  second-order  species  correlations  and  scalar  pdf's  in 
shear  layer  flows  and  should  provide  a  test  for  turbulence  modeling  cal¬ 
culations.  However,  as  discussed  later,  the  data  has  some  unexplained 
anomalies  that  reduce  its  value  for  model  verification. 

The  A.R.A.P.  multiequation  second-order  closure  code  for  shear  layer 
flows  has  been  used  to  calculate  the  flowfield  for  Helium  +  Argon  - 
Nitrogen  (constant  density)  and  Helium  -  Nitrogen  (density  ratio  *  7) 
planar  shear  layers.  The  velocity  ratio  for  both  flows  is  0.38.  The 
computations  are  started  at  the  exit  plane  of  the  splitter  plate  with 
smooth,  monotonic  initial  profiles  of  velocity  and  species  concentrations. 
The  initial  wall  boundary  layers  are  neglected  as  we  are  mainly  interested 
in  comparison  of  predictions  and  experimental  data  in  the  downstream 
fully  developed  region  of  the  flow.  A  triangular  profile  spot  of  turbu¬ 
lence  is  introduced  at  the  initial  station  with  u' u'  =  v' v'  *  w'w'  * 

2 | u ' v 1 |  .  The  peak  value  of  u' u'  *  O.Ol(Au)2  .  Previous  studies  have 


shown  that  with  reasonable  profiles,  downstream  results  are  quite  insen¬ 
sitive  to  the  width  and  amplitude  of  the  initial  turbulence  distribution. 


Calculations  have  been  carried  out  using  the  complete  "typical  eddy" 
pdf  model  including  density  fluctuations.  For  comparison  purposes,  cal¬ 
culations  were  also  made  using  a  pdf  "model"  in  which  all  third-  and 
higher-order  scalar  correlations  are  set  to  zero.  This  procedure  has 
been  termed  the  "second-order  approximation." 


Constant  Density  Shear  Layer  Studies 

For  nonreacting  constant  density  flows  the  transport  equations  for 
the  means  and  the  second-order  correlations  do  not  contain  any  third-order 
scalar  correlation  terms.  Therefore,  the  results  obtained  using  the 
"second-order  approximation"  and  the  "typical  eddy"  model  are  identical. 


Analysis  of  Data.  Brown  and  Roshko  (1974)  reported  measurements  of 
mean  velocity  profiles  in  air-air  mixing  layers  with  a  velocity  ratio  of 
0.38.  A  standard  measure  of  shear  layer  spread  rate  is  the  parameter  o0  . 
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y 2  and  y ^  are  the  profile  widths  corresponding  to  the  10%  and  90% 
velocity  profile  points  at  two  axial  positions  Xg  and  x^  .  u^  and 
Uj  are  the  freestream  velocities  in  the  two  streams.  We  have  calculated 
the  parameter  a0  from  the  Brown  and  Roshko  data  and  o0  =  9.7  .  This 


value  is  in  good  agreement  with  other  shear  layer  data  (Birch  and  Eggers, 
1973)  wherein  o0  ranges  between  9  and  11,  with  11  being  the  preferred 
value. 

Konrad  (1976)  reported  measurements  in  a  uniform  density  shear  layer 
composed  of  streams  of  He  +  Ar  and  Ng  in  the  same  apparatus  at  similar 
velocity  values,  and  o0  calculated  from  his  measurements  has  a  value  of 
only  7.5;  that  is,  the  shear  layer  is  spreading  too  rapidly.  This  dis¬ 
crepancy  in  the  spread  rate  between  Konrad  (1976)  and  Brown  and  Roshko 
(1974)  (and  other  data  in  the  literature)  has  to  be  resolved  before  one 
can  make  full  use  of  the  turbulence  correlations  and  pdf  measurements  of 
Konrad. 

Second-Order  Closure  Calculations.  The  predictions  for  the  mean 
velocity,  mean  species,  and  species  correlation  profiles  in  the  downstream 
similarity  region  are  compared  to  the  Konrad  (1976)  data  in  Figure  1. 

These  calculations  use  the  standard  A.R.A.P.  model  constants  that  have 
been  established  in  the  past  by  detailed  comparison  with  basic  incom¬ 
pressible  jet,  wake,  shear  layer,  and  flat  plate  boundary  layer  flows. 

The  multiequation  turbulent  reacting  £ode  (METREC)  uses  an  algebraic 
specification  of  the  turbulent  macroscale  A  . 
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Figure  1.  Mean  velocity,  mean  species,  and  species  correlation  profiles 
in  uniform  density  mixing  layer. 
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This  selection  of  the  length  scale  leads  to  a  prediction  of  aQ  =  9  . 

The  calculated  u  profile  is  somewhat  narrower  than  the  experimental 
profile  because  of  the  low  value  of  aQ  reported  by  Konrad. 

The  differences  in  the  mean  velocity  profile  for  this  flow  are  slight. 
However,  the  differences  in  the  mean  species  profiles  are  much  larger. 

To  match  the  experimental  data  of  Konrad,  it  is  necessary  to  increase  the 
scale  constant  C  to  0.7  and  to  reduce  the  turbulent  Schmidt  number  Scy  to 
0.3  from  its  standard  value  of  0.75  (Lewellen,  1977).  The  effect  of  these 
changes  on  the  calculations  is  shown  in  Figure  2. 

The  change  required  in  the  turbulent  Schmidt  number  is  quite  large 
and  is  not  supported  by  other  comparisons  with  data  for  constant  and  nearly 
constant  density  flows.  Thus,  Varma  et  al.  (1977)  have  demonstrated  good 
agreement  with  mean  temperature  axial  decay  rate  data  in  a  slightly  heated 
planar  jet  using  the  A.R.A.P.  models  and  $Cy  =  0.75  .  Lewellen  (1976) 
has  also  shown  good  agreement  with  measurements  in  a  heated  two-dimensional 
wake  using  a  scale  equation  formulation  and  Scy  =0.75  .  The  required 
value  of  Scy  =  0.3  seems  extremely  low  and  we  feel  additional  data  are 
required  to  justify  model  changes.  There  are  data  available  on  temperature 
profiles  in  a  thermal  mixing  layer  (Fiedler,  1974)  that  have  recently  come 
to  our  attention  and  this  should  provide  additional  information  on  the 
proper  value  of  the  scalar  transport  coefficients  in  a  turbulent  flow. 

The  computed  results  for  the  species  correlation  c'c'  in  Figure  2  show 
that  when  good  agreement  with  the  mean  quantities  is  secured,  the  fluctu¬ 
ation  level  is  also  in  reasonable  agreement  with  data. 
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y/x-x0 

Figure  2.  Mean  velocity,  mean  species,  and  species  correlation  profiles 
in  uniform  density  mixing  layer  with  increased  length  scale 
constant  and  decreased  turbulent  Schmidt  number. 
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This  selection  of  the  length  scale  leads  to  a  prediction  of  oQ  =  9  . 

The  calculated  u  profile  is  somewhat  narrower  than  the  experimental 
profile  because  of  the  low  value  of  aQ  reported  by  Konrad. 

The  differences  in  the  mean  velocity  profile  for  this  flow  are  slight. 
However,  the  differences  in  the  mean  species  profiles  are  much  larger. 

To  match  the  experimental  data  of  Konrad,  it  is  necessary  to  increase  the 
scale  constant  C  to  0.7  and  to  reduce  the  turbulent  Schmidt  number  Scy  to 
0.3  from  its  standard  value  of  0.75  {Lewellen,  1977).  The  effect  of  these 
changes  on  the  calculations  is  shown  in  Figure  2. 

The  change  required  in  the  turbulent  Schmidt  number  is  quite  large 
and  is  not  supported  by  other  comparisons  with  data  for  constant  and  nearly- 
constant  density  flows.  Thus,  Varma  et  al.  (1977)  have  demonstrated  good 
agreement  with  mean  temperature  axial  decay  rate  data  in  a  slightly  heated 
planar  jet  using  the  A.R.A.P.  models  and  Scy  =  0.75  .  Lewellen  (1976) 
has  also  shown  good  agreement  with  measurements  in  a  heated  two-dimensional 
wake  using  a  scale  equation  formulation  and  Scy  *  0.75  .  The  required 
value  of  Scy  =0.3  seems  extremely  low  and  we  feel  additional  data  are 
required  to  justify  model  changes.  There  are  data  available  on  temperature 
profiles  in  a  thermal  mixing  layer  (Fiedler,  1974)  that  have  recently  come 
to  our  attention  and  this  should  provide  additional  information  on  the 
proper  value  of  the  scalar  transport  coefficients  in  a  turbulent  flow. 

The  computed  results  for  the  species  correlation  c'c1  in  Figure  2  show 
that  when  good  agreement  with  the  mean  quantities  is  secured,  the  fluctu¬ 
ation  level  is  also  in  reasonable  agreement  with  data. 
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Variable  Density  Shear  Layer  Studies 


Analysis  of  Data.  There  is  again  a  significant  difference  between 
the  data  of  Brown  and  Roshko  (1974)  and  Konrad  (1976).  The  spread  rate 
parameter  o0  for  the  He  -  N2  shear  layer  as  calculated  from  Konrad  data 
has  a  value  of  7.4,  that  is,  virtually  the  same  o0  as  for  the  uniform 
density  shear  layer.  This  does  not  support  the  conclusion  of  increased 
spreading  rate  for  variable  density  flows  with  a  high  speed  Helium  stream 
configuration.  The  original  Brown-Roshko  data  predicts  oQ  =  8.8  , 
which  indicates  a  10-  15%  increase  in  the  spreading  rate  compared  to  the 
uniform  density  data. 

Second-Order  Closure  Calculations.  The  METREC  program  calculations 
have  been  carried  out  using  the  standard  model  constants  with  Scy  *  0.75 
and  a  length  scale  constant  of  0.6.  The  predictions  for  both  the  mean 
velocity  and  mean  species  profiles  are  significantly  narrower  than  the 
experimental  data.  The  calculated  spread  rate  a0  :  13  .  This  requires 
modifications  in  our  modeling  of  the  dynamics  of  the  turbulent  flow  and 
is  under  further  study.  The  important  variable  is  the  turbulent  macroscale. 
The  development  of  a  tensor  length  scale  equation  based  on  two-point  cor¬ 
relation  functions  is  underway  at  A.R.A.P.  as  part  of  other  research 
efforts,  and  the  same  procedure  should  be  useful  in  the  calculation  of  the 
turbulent  scales  in  variable-density  flows. 

We  have  also  carried  out  calculations  of  the  variable-density  shear 
layer  using  two-equation  turbulence  models  (ke2  model)  and  find  similar 
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disagreement  between  the  predicted  and  experimental  mean  profiles.  Parker 
and  Sirignano  (1978)  also  report  a  factor  of  two  error  in  the  predicted 
shear  layer  width  with  a  ke  model  compared  to  the  Brown-Roshko  data.  The 
ke  models  and  the  A.R.A.P.  second-order  closure  models  have  been  used 
successfully  for  computations  of  other  variable- density  jet  flows  (see 
below)  and  the  reasons  for  the  disagreement  with  the  shear  layer  data  have 
to  be  understood.  This,  of  course,  must  also  include  a  reexamination  of 
the  available  data  and  additional  experimental  measurements  in  variable- 
density  flows. 

The  results  for  the  species  correlations,  c'c'  ,  using  the  two  pdf 
models  are  compared  to  the  data  in  Figure  3.  The  error  in  the  mean  pro¬ 
files  has  been  taken  into  account  by  normalizing  the  results  by  the  width 
of  the  mean  species  profile.  The  results  for  c'c'  using  the  "typical 
eddy"  pdf  model  are  in  better  agreement  with  the  data  than  the  results 
obtained  by  using  the  "second-order  approximation"  procedure  that  neglects 
the  higher-order  scalar  correlations.  The  results  demonstrate  the  impor¬ 
tance  of  modeling  these  correlations  in  the  calculations. 

Variable- Density  Axisymmetric  Jet  Flow 

The  predictions  of  the  METREC  code  for  a  hydrogen  axisymmetric  jet 
flow  in  a  coflowing  airstream  are  compared  to  the  experimental  measure¬ 
ments  in  Figure  4.  The  experimental  data  is  by  Chriss  and  Paulk  (1972) 
and  is  Test  Case  10  of  the  NASA  Free  Turbulent  Shear  Flows  Conference  data 
set.  The  hydrogen  stream  has  a  velocity  of  1005  m/sec,  and  the  velocity 
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Figure  3.  Species  correlation  profiles  in  variable  density  mixing  layer. 


Figure  4.  Axial  decay  of  mean  velocity  and  mean  species  mass  fraction  in 
coflowing  axisymmetric  variable  density  jet  flowfield. 


ratio  of  hydrogen  to  air  stream  is  6.3.  The  density  ratio  at  the  nozzle 
exit  is  about  11.  The  calculations  used  a  turbulent  Schmidt  number  of 
0.75.  The  calculations  are  started  from  given  initial  mean  profiles  at 
x/D  z  3,  The  initial  turbulence  profiles  are  not  provided.  The  calcu¬ 
lations  are  started  with  a  spot  of  turbulence  at  the  initial  station. 

The  results  in  Figure  4  are  obtained  by  using  our  standard  procedure  of 
selecting  the  turbulence  level  u'u'  =  .01  Au2  at  the  initial  station. 

In  the  near  region  of  the  flow,  the  results  are  somewhat  sensitive  to  the 
initial  turbulence  level.  Figure  4  compares  the  predicted  centerline 
variation  of  the  mean  velocity  and  hydrogen  mass  fraction  with  the  data. 
The  agreement  is  quite  satisfactory.  For  this  flowfield  the  results 
obtained  with  the  use  of  the  "second-order  approximation"  are  virtually 
the  same  as  the  results  from  the  "typical  edC(y"  model  shown  in  Figure  4, 
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3.  STATISTICAL  CONSTRAINTS  ON  CORRELATIONS 
FOR  THREE  SPECIES  FLOWFIELDS 


Basic  statistical  principles  are  being  used  to  obtain  the  rigorous 
bounds  on  correlations  of  various  variables  in  a  variable  density  three 
species  flowfield.  The  procedures  being  used  are  similar  to  those  devel¬ 
oped  during  the  derivation  of  constraints  for  a  two  species  system  (Varma 
et  al.,  1978).  However,  the  three  species  system  is  far  more  complex  and 
at  present  we  have  only  succeeded  in  developing  some  of  the  bounds  for 
the  second-order  correlations. 

Consider  a  three  species  flow  with  constant  pressure  and  temperature. 
Let  a  ,  B  ,  and  y  represent  the  mass  fractions  of  the  species.  We 
are  mainly  interested  in  deriving  the  constraints  on  the  following  scalar 
moments:  a,  B  ,  ,  W  ,  aB  ,  op,  Bp,  and  p*"  .  Additional 

constraints  on  various  third-order  moments  will  also  be  useful. 

Given  a  Only 

We  first  assume  that  only  a  is  specified.  Simple  conservation 
relations  lead  to  the  following  inequalities. 


0  <  o  <  1 


0  <  B  <  1  -  a 


(1) 

(2) 
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Use  of  the  renormalization  theorem  and  Cauchy-Schwarz  inequality 


leads  to 


a  <  a2  <  a 


0  <  aB  <  a  -  a2 


0  <  W  <  3  -  a& 


Given  just  a  ,  Eq.  (3)  represents  the  tightest  bounds  on  a2  . 
Substituting  the  a2  bounds  in  Eqs.  (4)  and  (5)  leads  to 


0  <  a0  <  a(l  -  a) 


0  <  B2  <  1  -  a 


These  are  the  best  bounds.  The  extremum  values  correspond  to  the 


delta  function  pdf's  sketched  below. 


or  0 


0  5  1-5  1  0 


Gi  ven  a  and  6 

The  conservation  conditions  lead  to  the  following  expressions. 

0  <  a  <  1  -  I  (8) 

0  <1  <  1  •  a  (9) 

We  have  previously  derived  the  following  bounds  on  a2  and  l5"  , 


_ 2  _  _ 

a  <  a2  <  a 


(10) 


e2  <  e2  <  e 


(u) 


It  will  now  be  shown  that  pdf's  can  be  constructed  that  attain  these 
extremum  values.  Therefore,  the  above  represent  the  best  bounds  on  the 
second-order  moments. 

For  a  three  species  system,  the  allowed  phase  space  region  for  the 
pdf  is  represented  by  the  triangle  shown  on  the  following  page.  The 
specification  of  a  and  $  corresponds  to  defining  the  location  of  the 
centroid.  If  the  pdf  corresponds  to  all  the  mass  distribution  concen¬ 
trated  at  the  centroid,  the  lower  bounds  on  a5"  and  I5"  are  attained. 
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The  upper  bounds  are  attained  by  the  following  pdf's. 


=  0 


In  the  above  diagrams,  for  generality  we  have  considered  a  distri¬ 
bution  of  masses  along  the  lines  a  *  0  or  6*0,  respectively. 

These  mass  distributions  are  such  that  for  the  figure  on  the  left  the 
centroid  of  the  distribution  along  the  line  a  =  0  is  located  at  8*  , 
the  position  corresponding  to  a  straight  line  connecting  a  =  1  and  the 
given  centroid.  Similar  considerations  hold  for  the  figure  on  the  right. 
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The  bounds  on  ofi  are  more  interesting.  Appendix  A  details  the 
derivation  of  an  upper  bound  on  o0  .  Using  this  proof,  we  propose  the 
following  bounds  on  . 


o  <5e  <2-^4 

~  a  +  6 


(12) 


If  these  values  can  be  realized  by  physically  realistic  pdf's,  then 
these  are  the  tightest  bounds  on  the  correlation.  The  lower  bound  is 
attained  by  mass  distributions  that  are  restricted  to  the  lines  a  =  0 
and  6=0. 


a  /3  =  o 


The  upper  bound  is  realized  by  the  following  pdf. 


Therefore,  the  best  bounds  on  06  are  as  indicated  in  Eq.  (12). 
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Given  a  ,  B  ,  and  a2  (or  g2  or  aB  ) 

We  now  wish  to  determine  the  allowed  regions  for  the  (a^  ,  ]jT)  space 
and  for  the  (ct5”  ,  of)  space. 

The  independent  bounds  on  cJ  and  I5"  define  a  rectangle  in  the 
(a*  ,  W)  domain.  However,  for  certain  specified  values  of  the  second- 
order  moments,  there  are  regions  within  the  rectangle  that  must  be  ex¬ 
cluded. 


Consider  a2  =  a  :  To  attain  the  upper  bound  on  a*  ,  the  distri¬ 
bution  must  have  the  form  shown  on  the  following  page.  We  will  calculate 
the  smallest  value  of  W  that  is  allowed. 
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Cauchy-Schwarz  inequality  gives 


g2  > 

~  1  -  a 


21 


___ 

Thus,  fF  =  0  cannot  be  reached  for  maximal  a?  and  the  lower 
right-hand  corner  of  the  rectangle  is  excluded. 

Similarly,  by  reversing  the  roles  of  a  and  0  ,  for  0^  =  0  , 
the  minimum  allowed  will  be 


The  upper  left-hand  corner  of  the  rectangle  is  not  allowed. 


Consider  a2  =  a  :  We  wish  to  calculate  how  large  W  can  be  for 
the  lower  bound  of  a2  .  The  minimal  value  of  a2  is  attained  by  the 
distribution  sketched  below  on  the  left. 


s(a,  I -a) 


(5,0) 
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P(o,0)  =  P2(B)6(o-a) 


6  - 


02  = 


/' 

/' 


6P2(6)dB 


B2P2(B)dB 


consider  P2(B)  to  be  a  two  delta  function  distribution  with  the  delta 
functions  located  on  the  sides  of  the  triangle  (sketch  on  right). 

P2(e)  •  (l-w)6(B)  +  w6[B  -  (1-a)] 


B  *  w(l  -  a) 


B2  =  (1-a)2 

1-a 


W  *  B(l-a) 


(15) 


It  is  postulated  that  B7  cannot  be  larger  than  this  value  for 
other  mass  distributions  that  are  located  in  the  interior  region  of  the 
triangle.  This  has  not  yet  been  proven. 
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the  largest  allowed  a2  is 


_  _ 2 

Similarly,  for  $2  =  6  , 

a2  *  a(l  -  6)  (16) 

With  these  derivations,  the  excluded  sections  of  the  rectangle  are 
shown  below.  The  other  two  corners  of  the  rectangle  are  accessible  by 
pdf's  in  which  (i)  a  delta  function  is  located  at  the  centroid  (a,6)  and 
(ii)  three  delta  functions  are  located  at  the  three  vertices  of  the 
tri angle. 
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We  have  determined  the  shape  of  the  connecting  curves  between  the 
above  points  by  considering  a  series  of  two  delta  function  pdf's,  with 
the  delta  functions  located  on  the  sides  of  the  triangle.  The  line  join¬ 
ing  the  delta  function  locations  must  pass  through  the  specified  centroid 
(a,B)  •  The  allowed  region  has  been  calculated  for  the  special  case  a  * 

B  ■  0.25  ,  and  is  shown  in  Figure  5.  We  have  also  sketched  the  expected 
shape  of  the  allowed  region  for  a  =  0.3  ,  B  =  0.4  on  the  same  figure. 
There  are  some  indications  that  these  are  the  best  bounds  in  ( ) 
space,  but  we  have  not  yet  completed  the  proof  for  other  more  general 
mass  distributions  different  from  two  delta  functions  located  on  the 
sides  of  the  triangle; 

Similar  procedures  have  been  used  to  determine  the  allowed  region  in 
(a^.aB)  space.  The  independent  bounds  on  a*  and  a3  define  the  follow 
ing  rectangle  in  the  domain. 


Consider  aB  =  0:  We  wish  to  obtain  the  bounds  on  ar  .  For 
aB  *  0  ,  the  pdf  must  be  restricted  to  the  a  *  0  and  6*0  sides  of 
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statistical  bounds  on  (F  and  a7  for  specified 
a  and  B  in  a  three  species  system. 
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the  triangle.  We  take  a  pdf  composed  of  a  delta  function  located  at 
(a,0)  and  a  general  mass  distribution  on  the  line  a  =  0  . 

/3 

l 

£# 


P(a,e)  *  wdfcOPjU)  +  (l-w)6(e)6(a-a#) 
a  =  (1  -  w)a# 
a2  =  (1-  w)a2  *  a  •  a# 


(17) 


Consider  =  a  : 

P(a,e)  *  (l-a)P1(B)6(a)  +  o6(a-  1)6(3) 

5e*o  (18) 


The  right-hand  side  of  the  rectangle  is  not  allowed. 


Thus,  the  excluded  regions  of  the  rectangle  in  (a2,a6)  space  are 
shown  below 


From  Eq.  (4),  the  general  conservation  bound  on  aB  is 

fl  <  afi  <  a  -  b2 

For  given  a  ,  the  upper  bound  on  af  in  (ci^aB)  space  is  a 
straight  line  of  slope  -1 .  This  represents  a  line  joining  the  points 
and  (*»0)  •  The  upper  bound  is  attained  by  pdf  distri¬ 
butions  corresponding  to  two  delta  functions  located  on  the  sides  a  =  0 
and  y  ■  0  . 

The  shapes  of  the  connecting  curve  between  the  other  allowed  points 
can  again  be  numerically  calculated  with  the  simplifying  assumption  of 
pdf's  composed  of  two  delta  functions.  The  curve  connecting  the  two  top 
points  has  been  calculated  for  a  ■  B  ■  0.25  and  is  plotted  In  Figure  6. 
The  other  curves  still  have  to  be  calculated,  but  their  expected  behavior 
is  shown  in  the  figure. 
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o6  and  a5"  for  specified 
species  system. 


In  summary,  many  of  the  important  constraints  on  second-order  corre¬ 
lations  in  a  three  species  flow  have  been  derived,  and  will  prove  useful 
in  the  testing  of  three  species  pdf  models.  Constraints  on  some  correla¬ 
tions  involving  density  fluctuations  still  remain  to  be  developed  along 
with  the  completion  of  some  of  the  proofs  discussed  above. 
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4.  CONCLUSIONS 

Calculations  have  been  made  to  compare  the  predictions  of  a  second- 
order  closure  program  with  experimental  measurements  in  uniform  and 
variable  density  shear  layers  and  jet  flows.  The  studies  for  uniform 
density  shear  flows  show  good  agreement  for  the  mean  velocity  profiles 
and  shear  layer  spread  rate.  The  results  for  the  mean  species  profiles 
are  in  good  agreement  with  data  when  a  turbulent  Schmidt  number  of  0.3  is 
used.  This  is  significantly  lower  than  the  value  of  0.75  used  in  earlier 
studies,  that  are  in  good  agreement  with  data  for  planar  jets  and  wakes. 
Additional  data  on  scalar  profiles  in  shear  layers  will  be  useful  to 
resolve  this  discrepancy.  The  results  for  variable  density  shear  layer 
flows  indicate  the  need  for  some  improvement  in  the  modeling  of  the 
dynamics  of  turbulence.  The  models  appear  to  be  satisfactory  for  variable 
density  jet  flows.  The  modeling  of  higher-order,  scalar  correlations 
using  the  "typical  eddy"  pdf  model  results  in  significant  improvement  in 
the  prediction  of  second-order  species  correlations. 
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APPENDIX  A 


Consider  a  general  pdf  P(a,3)  of  the  form 

n 

P(ot,3)  =  ^w-jdCa-  ot.j)6(g-  3-j) 
i=l 


with  0  <  w-j  <  1  ,  0  <  <  1  ,  0  <  ^  <  1  ,  0  <  ai  +  3,  <  1  ,  and 


Then,  to  prove  the  inequality, 


a  3  <  =*-4 
a  +  $ 


(1) 


The  proof  is  by  induction. 
For  n  =  1  : 


P(a,3)  =  w5(a -  0^5(3-  3j) 
a  =  wa^ 


3  *  w3j 
a3  =  waj3j 


Al 


Substituting  in  Eq.  (1) 


al6l  £ 


alBl 

al+el 


or  aj  +  Bj  <  1  ,  which  is  true! 

We  assume  the  result  is  true  for  m  =  n-  1  . 

n-1 

Vi =  XWi 

i=l 


A2 


We  then  have 


n-1 


^](ai  +  Bi)wi 


1-1 

Then  we  wish  to  prove: 


Define 

n-1  n-1 

A  H  I*i*1  B  =  Z6JWj 

i=l  j*l 


a  *  cu  b  ■  c  =  wn 


Eq.  (2)  is  equivalent  to  showing 


AB 

A  +  B 


+  abc  < 


(A  +  ac)(B  +  be) 
A  +  ac  +  B  +  be 


(3) 


This  can  be  simplified  to 

(A  +  B)ab(A  +  B  +  ac  +  be)  <  A(Ab  +  abc)  +  B(Ba  +  abc)  (4) 

Define 

f(c)  =  (A  +  B)ab(A  +  B  +  ac  +  be)  -  A(Ab  +  abc)  -  B(Ba  +  abc) 

A  ,  B  ,  a  ,  and  b  are  independent  of  c  and  hence  f  is 
linear  in  c  . 

To  prove  inequality  (4)  is  equivalent  to  showing  that 
f(0)  <  0  and  f'(0)  <  0  . 

f (0)  *  (A  +  B)2ab  -  A2b  -  B2a 


Lenina  1 


f(0)  <  0 


(A  +  B)2ab  <  A2b  +  B2a 


(5) 


A4 


Proof  of  Lemma 


If  a 
Assume 

Let  x 


or 

or 

Now  to 

or 


or  b  is  zero,  the  conclusion  is  trivial, 
a  j*  0  ,  b  >*  0  -*•  a<l  ,  b<l 


•  I  V  .1 

a  ’  y  b 


a  <  1  -  x  =  1  +  ex 

for  some 

Cj  >  0 

b<l  y=l  +  e2 

for  some 

Eg  >  0 

a*b<  1  *  i  ♦  7<1 


1  4  <  1 


1  +  Ej  1  +  Cg  - 


1  i  C1 '  e2 


prove, 


(A  +  B)2  <  A2(l  +  Cj)  +  B2(l  +  e2) 


2AB  <  A*e1  +  B2c2 


A5 


consider. 


(A  -  B  t/e^)2  >  0 


A2ex  +  B2e2  >.  2AB  /e jCg 

>  2AB  as  /ejE^  >  1 

This  proves  the  lemma. 

Lemma  2 

f'(0)  <  0 

f'(0)  *  (A  +  B)ab(a  +  b)  -  Aab  -  Bab 
f'(0)  *  ab(A  +  B)(a  +  b  -  1) 
a  +  b  -  1  <  0 


f'(0)<0 

The  proof  is  complete. 

of  <  ~ 
a+  B 
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